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Abstract 

We analyze the potential sensitivity of a search for e — )• r conversion at a proposed electron- 
ion collider (EIC) facility. To that end, we calculate the cross sections for e — )• r events in a 
leptoquark framework assuming that the leptoquark masses are on the order of several hundred 
GeV or more. Given present limits on leptoquarks from direct searches at HERA and rare decay 
processes, an EIC sensitive to 0.1 /6 e — )• r cross sections could probe previously unexplored regions 
of parameter space for these lepton flavor violating events (assuming 90 GeV center-of-mass energy 
and 10 fh~^ integrated luminosity). Depending on the species of leptoquark and flavor structure 
of the couplings, an EIC search could surpass the HERA and rare process sensitivity to e — r 
conversion amplitudes by as much as an order of magnitude or more. We also derive updated limits 
on quark flavor-diagonal LFV leptoquark interactions using the most recent BaBar r — t- 67 search. 
We find that limits from an EIC e — t- r search could be competitive with the most recent r — t- 67 
limit for a subset of the quark flavor-diagonal leptoquark couplings. Using an SU(5) GUT model 
in which leptoquark couplings are constrained by the neutrino masses and mixing, we illustrate 
how observable leptoquark-induced e — t- r conversion can be consistent with stringent LFV limits 
imposed by — )■ 67 and /i — )■ e conversion searches. 
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I. INTRODUCTION 



The possibility that flavor may not be conserved by charged leptons at an observable 
level continues to be a topic of substantial interest in particle physics. Although neutrino 
oscillations imply that charged lepton flavor violating processes such as /i — 67 are allowed 
in the Standard Model (SM), their rates are extraordinarily small — and therefore unob- 
servable — because they are suppressed by the masses of the neutrinos. Many models of 
physics beyond the SM predict rates for charged lepton flavor violation (LFV) processes 
that are both larger than the SM predictions and within reach of present or future experi- 
ments. Hence, LFV is considered an important probe of new physics models such as SUSY 
SO{10) grand unification theories (GUTs) (for a concise review of several models, see (H), 
a supersymmetric SM with right-handed neutrinos j2|, left-right symmetric models and 
the Randall- Sundrum model (see for example 0] and the references therein). LFV is also 
present in models containing leptoquarks with trans-generational couplings to leptons and 
quarks. Various types of leptoquarks are predicted to exist, for example, in the Pati-Salam 
SU{4) color model [sl and non-supersymmetric GUTs utilizing SU{5) jof and 5*0(10) 0| 
symmetry groups. 

Extensive searches for charged lepton flavor violation between the first and second lepton 
generations (which we will refer to as "LFV(1,2)" for brevity) have placed stringent experi- 
mental limits on processes such as /i — ?■ e conversion (from the SINDRUM II collaboration 
Q), /i ^ 67 (from the MEGA collaboration [loj), and /i ^ 3e (from the SINDRUM 
collaboration ll|)lll LFV in the first and third generations ("LFV(1,3)") is also possible, 



although current experimental limits on r — )■ 67 (from the BaBar collaboration 13|), and 



r — )■ 3e (from the BELLE collaboration IJ]) are several orders of magnitude weaker than 
their LFV(1,2) counterparts. The gap between the LFV(1,2) and LFV(1,3) limits will con- 
tinue to widen with the next generation of experiments searching for LFV (1,2), such as 



Fermilab's proposed Mu2e experiment [15j and the MEG experiment at PSI [16|. 

In this study, we examine the prospects for LFV(1,3) searches at a prospective electron- ion 
collider. Given the current status of experimental limits just described, one may ask whether 
incremental (~ two orders of magnitude) improvements in LFV(1,3) sensitivities are still 
useful as a probe of new physics. The answer is certainly yes if there exist models in which 
processes such as r — > 67 and e ^ t are enhanced by several orders of magnitude relative 
to u — >■ 67 and e — >■ (or — )■ e) conversion. One example of such a model is discussed in 



17[ , in which the authors use a particular parameterization of the minimal supersymmetric 
seesaw model to find regions of parameter space which suppress the /i — 67 branching 
fraction and enlarge the r — 67 branching fraction. Another model in which LFV(1,3) can 



be enhanced is the SU (5) GUT with leptoquarks introduced in [18| and further studied in 
19| . In this model, the leptoquark couplings are tied to the neutrino masses and mixing 
matrix, with the result that in the quasi-degenerate neutrino mass regime the cross section 
for e — )■ r can be relatively large despite strong constraints on /i — )■ e (we will describe this 
analysis in a later section of the paper). In light of these scenarios, the theoretical motivation 
for our analysis is clearly present. 

In what follows, we analyze the prospective physics reach of a search for LFV(1,3) at 
the proposed electron-ion collider (EIC), a high luminosity accelerator facility with the 
primary goal of exploring several open questions in QCD, such as the properties of the 

^ We cite the experiments with the strongest hmits. Other searches for these processes can be found in the 



Particle Data Group hstings 12 1 
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gluon distribution in nucleons and high-density quark-gluon matter [20j. In the eponymous 
electron-ion colhsions at the EIC, the LFV(1,3) process of interest is e — )■ r conversion. The 
e — >■ r process can occur in leptoquark models via tree-level interactions, and so stronger 
signals may be expected relative to other models of LFV(1,3) that induce lepton flavor 
violation through loop effects (e.g., doubly-charged Higgs loops in the left-right symmetric 
model). Therefore, in our initial analysis of LFV(1,3) at the EIC presented in this paper, 
we consider searches for leptoquark-induced e — )■ r events. For completeness, we adopt the 
general BuchmuUer-Ruckl-Wyler leptoquark parameterization, described in section [Tll 

As we will discuss in section IIIH the ZEUS and HI experiments at HERA searched for 
e — )■ r leptoquark events and placed upper limits on the ratio of the leptoquark-quark-lepton 
couplings divided by the squared leptoquark mass. In addition to the direct searches for 
leptoquarks at HERA, experimental limits on r — )■ 67 and other rare decays (e.g., r — )■ vre, 
T — > 3e, and decays of K and B mesons) also place constraints on the same couplings-over- 
mass ratios [2l|. At the time the HERA results were published, the limits from r — )■ 67 were 
weaker than limits from other rare processes and were not relevant for the analyses by the 
ZEUS and HI collaborations. In 2010, a stronger limit on the r — )■ 67 decay was published 
by the BaBar collaboration [l3[, and so in section [TV] we update the leptoquark limits using 
this most recent experimental result. 

Given an EIC operating at 90 GeV center-of-mass energy and with 10 fh^^ of integrated 
luminosity, we will show in section |V] that the EIC could set limits on the leptoquark 
coupling-over-mass ratios that surpass the current best limits from the HERA experiments 
by as much as nearly two orders of magnitude. We will also show that the EIC could compete 
with or surpass the updated leptoquark limits from r — )■ 67 for a subset of the quark flavor- 
diagonal leptoquark couplings. Given this potential physics reach, it is interesting to ask 
whether there exist theoretically well-motivated scenarios in which e — )■ r conversion could 
be observed at the EIC despite the much stronger limits on LEV from /i — )■ 67 and /i — ?■ e . 
To that end, in section IVIl we discuss the interesting connection between one of the BRW 



leptoquarks and the model presented in [18|. The leptoquark of interest is unconstrained 
by yU — )■ 67 and r — )■ 67, and using the results of the analysis in 3] we show that this 
leptoquark can yield e — )■ r cross sections within reach of the EIC and still be compatible 
with the limits on neutrino masses and mixing angles and /i — e conversion. In section IVIIl 
we will summarize our results and comment on a few considerations relevant for undertaking 
an experimental search for e — > r at the EIC. We will briefly make the point that searches 
for leptoquarks in e — )■ r events at the EIC are complementary to leptoquark searches at 
hadron colliders like the Tevatron and LHC. 



II. LEPTOQUARK FRAMEWORK 



We use the Buchmuller-Ruckl-Wyler (BRW) leptoquark parameterization 22] as the 
framework for our LFV(1,3) analysis. The BRW parameterization catalogs all possible 
renormalizable and S'f/(3)c x SU{2)l x U{1)y invariant interactions between (scalar or 
vector) leptoquarks and SM fermions. These interactions are given by the Lagrangian in 
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eq. dZT 

'■LQ = i 



C-LQ — C,F=0 + C,\F\=2 



+ h^dRl.enV,''^ + h^UR^^CRV^''^ + /if g.^TM^^LV^i''^ + h-c (2.1) 

In eq. (12. ip . and are the SU{2) doublet quarks and leptons, ur, (Ir, cr are the SU{2) 
singlet quarks and charged lepton, e is the SU{2) antisymmetric tensor (ei2 = —621 = +1), 
f = (ti, T2, Ts) are the Pauli matrices, and the charge conjugated fermion is defined as iJj'^ = 

-F T 

Cip = ^7270'^ in the Dirac basis for the 7 matrices. Color, SU{2), and fiavor (generation) 
indices have been suppressed. The leptoquarks are characterized by their fermion number, 
their spin, the chirality of their coupling to leptons, and their gauge group quantum numbers. 
The leptoquarks carry fermion number F = 3B + L equal to or ±2. We follow the notation 
used in the recent hterature where spin-0 leptoquarks are S and spin-1 are V, the subscript 
indicates the SU{2) quantum number (0 for a singlet, 1/2 for a doublet, 1 for a triplet), 
the superscript L, R indicates the chirality of the lepton coupling to the leptoquark, and a 
tilde (~) is used to distinguish between leptoquarks which have different hypercharges but 
are otherwise identical. The dimensionless coupling constants g and h (which we assume to 
be real) carry the same lepton chirality and SU{2) labels as their associated leptoquarks. 
Lepton fiavor violation can arise if the couplings — which are matrices in fiavor space — 
have non-zero off-diagonal elements. 

We will also require the interactions between the BRW leptoquarks and the photon. The 
photon interactions arise from the Lagrangian kinetic terms with SU(2)l x U{1)y covariant 



derivatives acting on the leptoquark fields (23[: 



^ir^Z^ = (D.SnD^S) , (2.2) 



rirlrJ = -I (D^V, - D,V,)^ (D'^V - D'^Vn . (2.3) 



^kinetic ' 2 



D, = d^, + igf-W, + ig'-B^ , (2.4) 



The covariant derivative is given by 

.Y 
2 

where the T"" are the generator matrices for the SU{2) representation occupied by the 
leptoquarks (singlet, doubletH, or triplet). The photon interaction for a scalar leptoquark is 
given by 

4qT = '^Q^Q ii^^') ^-S' (9,S)] , (2.5) 

where Qlq is the electric charge of the leptoquark. 

For the vector leptoquarks, interactions with the photon depend on the nature of these 
massive vector particles, i.e., whether or not the leptoquarks are gauge bosons of some 



^ Note that the doublets must be in the 2 representation given the form of the Lagrangian in eq. (|2.ip . 
E.g., exphcitly writing the SU{2) indices, URi^LiSy^ ■ shows that the i — 2 component of the leptoquark 
multiplet couples to the electron and must have the opposite eigenvalue to be SU{2) invariant. 
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beyond-the-SM symmetry group. In addition to the interaction arising from eq. (12. 3p . there 
can exist an anomalous magnetic moment couphng of the leptoquark to the photon, so the 
full interaction Lagrangian is 



^{vector) 



where the leptoquark field strength tensor V^*^ is given by 



(2.6) 



(2.7) 



and F'^'^ is the usual photon field strength tensor. If the leptoquarks are gauge bosons (as 
in the case of some SU{5) GUTs, e.g.), then n = and the resulting photon interaction is a 
three-gauge-boson vertex, the result of spontaneous symmetry breaking of the higher gauge 
group containing both the leptoquarks and the photons to U (1) em- (Also, if the leptoquarks 
are gauge bosons, eq. (12. 3 p is replaced by the appropriate kinetic term for the gauge bosons 
of the larger symmetry group.) This question of the gauge nature of the vector leptoquarks 
will have further implications for our analysis, particularly in the calculation of the r — )■ 67 
limits (see section lIVp . Finally, the electric charges of the scalar and vector leptoquarks 
which appear in the photon interaction terms are easily determined from eq. (12. ip (also, see 



Table 1 in |23|). 



III. CROSS SECTION CALCULATIONS FOR e ^ r 

Electron to tau conversion in an e~p deep inelastic scattering process is the LFV(1,3) sig- 
nal at the EIC which we consider in our analysis. In the BRW leptoquark parameterization, 
such a process occurs via tree level partonic interactions. In e~p collisions, F = type lepto- 
quarks couple to antiquarks in the s-channel and quarks in the u-channel, while |-F| = 2 type 
leptoquarks couple to quarks in the s-channel and antiquarks in the w-channel (see fig. [1]). 
If the leptoquark mass is much larger than the center of mass energy, Mlq ^ y/s, the 
momentum dependence of the leptoquark propagator can be neglected, effectively shrinking 
the partonic interaction to a four-fermion vertex. The cross section then depends only on 
the ratio of the leptoquark couplings divided by the leptoquark mass. The total inclusive 
cross section for e~ + p ^ t" + X with a single intermediate leptoquark is given (in the 
limit of massless quarks and leptons) by HT 



O'F=0 



E 

a,l3 



(T\F\=2 = 



S 

32^ 

s 

32^ 





to 












to 











dxdy xq^ {x, xs) f (y) + J dxdy xq^ (x, -u) g (y) 
dxdy xqa {x, xs) f {y) + / dxdy xq^ (x, -u) g (y) 



(3.1) 



The functions / and g are defined in eq. (13. 2p . 



fiy) 



1/2 (scalar) 
2(1 — y)'^ (vector) 



9{y) 



/I (scalar) 
2 (vector) 



(3.2) 
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s-channel u-channel s-channel n-channel 

F = 1^1 = 2 



FIG. 1. Feynman diagrams for e — )■ r scattering processes via leptoquarks which depend on the 
parameter AiaAa^/M^g. The partonic cross section is convoluted with the pdf of the initial state 
(anti)quark of each diagram. See eq. (|3.1|) . 

The parton distribution functions for the quarks and antiquarks are q (x, Q'^) and q (x, Q"^), 
respectively, evaluated at momentum fraction x and energy scale Q^. Also, u = xs [y — 1) 
and both x and y are integrated from to 1. The leptoquark couplings Aiq and As^ are 
the couplings g and h which appear in the Lagrangian of eq. fl2.1l) (additional factors of — 1 
and/or a/2 may multiplythese couplings, depending on the leptoquark SU{2) representation 
— see, e.g.. Table 2 of |22| and Table 1 of 23|). The subscripts on the couplings A are gener- 
ation indices: 1 and 3 for the electron and tau, and a and (3 for the quarks/ant iquarksjfl We 
refer to ratios with a = /3 as "quark flavor-diagonal" and those with a 7^ /3 as "quark flavor- 
off-diagonal" . The ZEUS and HI collaborations placed upper limits (at 95% confidence 
level) on the ratio XiaXsp/M^Q for each type of BRW leptoquark and for all combinations 
of a and /3 except in cases where the top quark was the only third-generation quark coupled 
to the leptoquark (23-27|. To obtain these limits, several assumptions were made: only 
one type of leptoquark dominated the cross section, the leptoquark coupled only to left- or 
right-handed leptons but not botli^, and leptoquarks in SU{2) multiplets are degenerate in 
mass. We make these assumptions in our analysis as well. 

To determine the sensitivity of an EIC search for LFV(1,3) in e — )■ r processes, we 
calculate an upper bound on the cross sections for the various leptoquarks using eq. (13. ip 
and the most stringent limits on Aiq,A3/3/M|q from the ZEUS or HI collaborations (or those 
rare process limits cited by the ZEUS and HI collaborations). We use the MSTW 2008 
NLO set for the quark and antiquark proton p.d.f.sH From eq. (13. ip . there is a one-to- 
one correspondence between the partonic sub-process cross section and the leptoquark ratio 
AiaAs^/Mlg. Given a number for the sub-process cross section, we calculate the leptoquark 
ratio and then scale (i.e., divide) the leptoquark ratio by the HERA/rare process limit. We 
define this scaled leptoquark ratio as the variable z. Thus, for a given cross section there 
is a unique value of z. In other words, z is the fractional reduction in the leptoquark ratio 
relative to the HERA/rare process limit. Results of these calculations will be presented in 
section IVl after we discuss limits from r — 67. 

^ Note that a is not always the initial state quark/ antiquark; sec fig. [H 

^ This assumption was already made in writing the Lagrangian in eq. (j2.ip . Leptoquarks with identical 
quantum numbers, e.g. Sq and S'(f , could have identical couplings to left- and right-handed leptons: 
9o = 9o 

leptons were not differentiated. 
^ http: / /projects. hepforge.org/ mstwpdf/ 



L _ „R. the original BRW parameterization [22| , leptoquarks coupling to both left- and right-handed 
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FIG. 2. Leptoquark loops contribute to r — ?■ 67 decay. 



IV. UPDATING T ^ 67 LIMITS 



As mentioned in the introduction, the BaBar collaboration has published a stronger limit 
(l3| on the branching fraction for r — )■ 67 since the time of the ZEUS and HI analyses 
(which did not make use of the weaker contemporaneous r — > 67 limit). The experimental 
bound on r — )■ 67 implies additional constraints on the ratios Ai^As/j/Mlg independent of 
the limits from HERA and other rare processes such as r — j- 3e, r — ire, and B and K 
decays. However, the r — )■ 67 bound only applies to those leptoquark ratios where a = (3 
( "quark flavor-diagonal" ) since the r — > 67 process proceeds via a loop with a single quark 
and a leptoquark, shown in fig. [2J In this section, we will use the recent limit from BaBar 
to calculate new leptoquark limits from r — >■ 67. 

In general, the amplitude for a r — >■ 67* process can be written 0, 28 



M 



--ee*'u, ip') [{q^^, - q^) {A^Pl + A^Pr) 

+imrq''a^a {A^Pl + A^Pr)] (p) 



(4.1) 



For a real photon (g^ = 0), only the magnetic moment term containing the coefficients A2 
and A2 will contribute to |A^|^; the branching ratio for the r — 67 process is then 

r (r~ e~7) ASn^aEM 



Br (r — 67) 



r (r^ —7- e'lygUT-) 



G2 



A 



L|2 



A. 



i?,|2 



(4.2) 



In eq. (14. 2p . is the Fermi constant obtained from muon decay. 

We calculate the contributions to A^ and A2 for each scalar leptoquark (vector lep- 
toquarks will be discussed later) by computing the amplitude for the diagrams in fig. [2] 
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according to the Feynman rules obtained from the leptoquark Lagrangian (eq. fl2.ip ) and 
from the photon interactions in eq. (12. 5p . In the hmit of zero electron mass, only the first 
two diagrams of fig. [2] are needed; however, all four diagrams are necessary to cancel the 
loop divergences and thus avoid introducing a tree level counterterm for a r, e, 7 vertex. To 
evaluate the loop integrals, we approximate m^/Mf^Q ^ and expand in powers of m^/MlQ. 
With this procedure, our results are valid for all quarks running in the loop, including top 
quarks. For all of the scalar leptoquarks, the expressions for A2 and A2 have the same basic 
structure: the A2 coefficient with the same left- or right-handed label as the leptoquark is 
zero, while the other coefficient is given by 




2 

11 + 6 In ^ 




In eq. f l4.3p . A''^ = 3 is the number of colors and Qg (Qlq) is the sum of the electric charges 
of all the quarks (leptoquarks) appearing in the loop — note that leptoquarks in SU(2) 
doublets or triplets can couple to both up- and down-type quarks, with different electric 
charges for each member of the multiplet. The quark charge contribution is from diagram (a) 
of fig. El while the leptoquark charge contribution comes from diagram (b). The leptoquark 
charges are determined by their SU{2) eigenvalue and hypercharge, Q = T3 + Y/2, and 
can be inferred from the Lagrangian in eq. (12.11) . Complete results for the individual scalar 
leptoquarks are given in appendix \M We note that our results for the scalar leptoquarks, 
eq. (143 p . agree with the results in ]2l| and j29|FI 

The uncertain gauge nature of the vector leptoquarks presents a difficulty in calculating 
the r — )■ 67 loop amplitudes. In particular, the interaction between the photon and lepto- 
quark depends on the unknown parameter k. Moreover, if the vector leptoquarks are gauge 
bosons (k = 0), their propagator can be written (in the 't Hooft-Feynman gauge) as 



P - MIq + le 



(4.4) 



but if the vector leptoquarks are not gauge bosons, then their propagator is 



The second term in eq. (14. 5 p introduces extra divergences in the loop graphs of fig. |2] that 
do not cancel like the divergences of the scalar leptoquark diagrams. The authors of 21 
argued that the second term in the propagator of eq. (14.51) can be neglected for the purpose 
of extracting upper limits on the leptoquark coupling-to-mass ratios: since the additional 
divergences introduced by this term go like powers of a large (model-dependent) cutoff scale, 
ignoring the extra contribution to the propagator results in conservative (i.e., weaker) upper 
limits. 

Here, we adopt a different perspective. The manifest non-renormalizability of the massive 
vector leptoquark theory leads to several options. It is possible, for example, that the vector 



^ The authors of ull assumed massless quarks and therefore do not have a term proportional to m^/Mlg. 
The author of [29| assumes that the couphngs A arc unitary, so the first term in eq. (|4.3p which is 
independent of the quark mass vanishes when summing over the quark generations. 
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leptoquarks are gauge bosons of a larger gauge group and that they receive their masses 
through spontaneous breaking of the gauge symmetry down to that of the Standard ModeL 
In this context, one would presumably need to include loop contributions involving other 
massive degrees of freedom associated with the extended gauge group, and it is not clear 
without specifying a model whether these contributions will add to, or cancel against, the 
leptoquark loops. Alternatively, one might consider the leptoquarks (not necessarily as 
gauge bosons) as part of a low-energy effective theory valid below a scale A ~ Mlq. In 
this case, there would exist a set of higher dimension LFV operators with a priori unknown 
coefficients determined by physics above the scale A. In this situation, it is possible to 
derive "naturalness" bounds on the leptoquark couplings by requiring that the individual 
contributions from cut-off dependent loop amplitudes and higher dimension operators be 
no larger than the r — 67 bounds. Depending on which option one chooses, one may 
derive different limits on the coupling-to-mass ratios. In short, it is not possible to derive 
theoretically robust constraints using the particle content of eq. fl2.ip alone. Therefore, we 
do not attempt to infer r — 67 limits for vector leptoquarks and restrict our attention to 
the scalar lepto quar ks. 

Experiments jl2l . 13| have given 



^ ^"'^^'^-^ - 0.1785 , < 3.3 X 10-« (4.6) 



r (r total) ' r (r total) 

so we take Br (r — 67) < 1.85 x 10"'^. From eq. (14. 2 p and eq. (14. 3p . it is evident that this 
experimental upper bound on the branching ratio for the r — )■ 67 decay places an upper 
limit only on a sum over generations of the leptoquark ratios )^ia^za/ ^Iq'-i furthermore, the 
terms proportional to the quark mass contain additional powers of the unknown leptoquark 
mass in the denominator. So, we will first neglect the terms proportional to m^/Mlg since 
we assume the leptoquarks have large masses, Mlq ^ C (100 GeV). The quark mass 
is expected to make a significant contribution to the A2 coefficients only if the quark in 
question is a top quark, and we will comment on the effect of including the top quark mass 
at the end. Then, to extract an upper limit on the ratio for a single generation, we consider 
two different options. The first option is to assume that each individual quark generation 
saturates the r — )■ 67 limit. This assumption gives a weaker (larger) upper limit on the ratio 
for each generation a. The second, more "democratic", option assumes that all generations 
contribute equally to the r — >■ 67 branching ratio, giving a stronger (smaller) upper limit 
on the leptoquark ratio for each generation. In both options, the limit is equivalent for all 
generations. These two options determine a range of values for the upper bound on the 
leptoquark ratios XiaXsa/Mlq arising from the experimental r — )■ 67 limit. Table |T] gives an 

example of these limits for the leptoquark Sq. 

Including the top quark mass will result in an increase in the calculated upper bound 
on Aiq,A3q,/M|q; the increase is mostly independent of the leptoquark type (since all the 
leptoquarks coupling to the top quark have charges that are O (1)) but dependent on the 
leptoquark mass itself and which of the two options is used to calculate the upper bound 
on XiaXsa/Mlq. We show in table Hi] the percent increases in the leptoquark ratio upper 
bounds under various assumptions. By neglecting the top quark mass, we get stronger upper 
limits on the leptoquark ratios from r — ?■ 67; i.e., the ratios Aiq,A3q,/M|q are smaller when 
the top quark mass is neglected. This is preferable for making comparisons with the ability 
of the EIC to probe these leptoquark ratios since smaller limits from r — )■ 67 will force more 
conservative (larger) estimates of the integrated luminosity necessary for the EIC to surpass 
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Option 1 


Option 2 




0.481 


0.160 TeV~'^ 


a = 1, z = 


1.2 


0.40 


a = 2, z = 


0.074 


0.025 


a = 3, z = 


0.032 


0.011 



TABLE I. r — 67 upper limits on AiqAsq/M^q for the leptoquark Sq. Option 1 assumes each 
generation individually saturates the r — )• 67 limit, while Option 2 assumes all three generations 
have equal contributions. We compute the scaled leptoquark ratio z by dividing the r — )• 67 limit 
on XiaXsa/M^Q by the limit from HERA and other rare processes, z > 1 implies that the bound 
from HERA/rare processes is already stronger than the new limit from r — )• ej. These limits 
appear as the vertical dashed lines in fig. [3] and fig. 21 





Mlq = bOOGeV 


Mlq = ISOOGey 


Option 1 


+25% 


+10% 


Option 2 


+7% 


+3% 



TABLE II. Estimates of the effect of including the top quark mass in calculating the upper limits 
on Xia^sa/Mlg from r — )■ ej. Including the top quark mass weakens (i.e., increases by the 
percentage indicated) the upper bounds on the leptoquark ratios Aiq,A3q,/M|q calculated from the 
T — )■ 67 experimental bound. The top quark mass effect depends on the leptoquark mass and the 
details of how the leptoquark ratio AiqA3q,/M|q is calculated (described in the text as Option 1, 
in which a single quark generation saturates the upper bound, and Option 2, in which all three 
quark generations contribute equally to the upper bound). 

the r — >■ 67 limits. In the next section, we combine the r — )• 67 limits with the cross section 
calculations discussed in the previous section and present our results. 



V. RESULTS AND DISCUSSION 

We now present our numerical results and discuss their implications for the EIC. In 
principle, with 10 fb~^ of integrated luminosity, the EIC could probe previously unexplored 
regions of parameter space for leptoquark-induced e — )■ r events with a cross section of 0.1 
fb or greater. For the specific case of e —> r events involving leptoquarks, we make use 
of eq. fl3.ip to plot the partonic sub-process cross sections for individual BRW leptoquarks 
as a function of z, the leptoquark ratio XiaXsp/Mlq scaled by the corresponding HERA 
limit. A center-of-mass energy of a/s = 90 GeV (roughly corresponding to, e.g., a 10 GeV 
electron beam and 200 GeV proton beam) is assumed. We separate the contributions to the 
total inclusive cross section a {e~p — )■ t~X) from all combinations of initial and final state 
quark/antiquark generations; the separate contributions are proportional to AiqA3/3/M|q 
where a and (3 are the quark generation numbers, as discussed above. A value of 2; = 1 
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FIG. 3. Cross section for the scalar leptoquark Sq plotted according to eq. (|3.ip as a function of z, 
defined in the text as the leptoquark ratio Xia^3p/M^Q divided by the HERA/rare process limit on 
that ratio. This plot shows the cross section with second generation initial and final state quarks, 
(a, /3) = (2,2). Upper bounds on A12A32/MIQ are calculated from the r — )• 67 limit according to 
the two options discussed in the text; these upper bounds are indicated by the vertical dashed lines 
joined by the horizontal arrow. The upper limit on A12A32/MIQ from r — )• 67 may be as large as 
the vertical line on the right or as small as the vertical line on the left. 

corresponds to the maximum cross section allowed by the limits from HERA El For scalar 
leptoquark quark flavor-diagonal cases where a = /3, the r — )■ 67 upper limits on XiaXsa/Mlg 
derived according to the two options discussed in section |IV] are also applicable; these bounds 
are also scaled by the corresponding HERA limit to obtain values for the variable z. 

In fig. [31 we show for the particular leptoquark Sq the cross section calculated according 
to eq. (13. ip as a function of the leptoquark ratio AicAs^/Mlg (scaled by the corresponding 
HERA limit to give the variable z). We have chosen (a, /3) = (2, 2), corresponding to second 
generation quarks/antiquarks in the initial and final states (see fig. [T]). Since a = /3 is 
"quark flavor-diagonal", limits from r — >■ 67 apply: the vertical dashed line on the right 
corresponds to the upper bound on A12A32/MIQ calculated according to Option 1 and the 
vertical dashed line on the left corresponds to Option 2, as discussed in section IIVI The 
upper bound on the r — j- 67 branching ratio corresponds to an upper limit on X12X32/ MIq 
that may lie anywhere between the two vertical dashed lines (this region is indicated by 
the horizontal arrow). The r — > 67 bounds on the leptoquark ratios have also been scaled 
by the HERA limit. The horizontal dashed line indicates a 0.1 fb cross section, which is 
roughly the level of sensitivity the EIC could achieve with 10 fb~^ integrated luminosity. 
If no e — r events were observed, the EIC could set a new limit on the ratio X12X32/ M^q 
given by the intersection of the diagonal line with the horizontal dashed line. This occurs 



^ We use the most restrictive value of XiaXsp/M^Q from the most recent ZEUS [25[ and HI 27[ analyses 
and the rare process results cited therein. In fact, for all e — >■ r leptoquark processes, the ZEUS limits are 
stronger than the HI results, although this is not the case for the rare process limits. The rare process 
limits used in the HERA analyses did not make use of the relatively weak bound on r — >■ 67 which existed 
at the time. 

11 




FIG. 4. Cross sections for various scalar leptoquarks plotted according to eq. ()3.ip as a function 
of the scaled leptoquark ratio AiqAs/j/MIq. These plots show the quark flavor-diagonal cross 
sections with (a,/3) = (1,1) ,(2,2) ,(3,3). The range of the leptoquark ratios Aiq,A3q,/M|q that 
satisfy the r — )• 67 limit are indicated by the vertical dashed lines and horizontal arrows. Cross 
sections and r — )• 67 limits for third generation quarks are not shown if the leptoquark couples 
exclusively to the top quark; no limit from the HERA analyses exists for these cases. 




FIG. 5. Cross sections for various scalar leptoquarks plotted according to eq. (13. ip as a function 
of z. These plots show the quark flavor-off-diagonal cross sections with a ^ (3. Cross sections for 
third generation quarks are not shown if the leptoquark couples exclusively to the top quark; no 
limit from the HERA analyses exists for these cases. 



12 



for z ^ 0.1; hence, the new hmit would be one order of magnitude smaller than the current 
limit. The new limit on A12A32/MIQ would surpass the weaker r — ?■ 67 upper bound from 
Option 1, but not the stronger limit of Option 2. 

In fig. m and fig. we show the cross sections calculated for four different scalar lep- 
toquarks, 5'^/2' '^1/2' '^(f ' ^'^'^ "^(f ' ^ function of the scaled leptoquark ratio, z. We have 
chosen to display results for these four leptoquarks because they demonstrate many gen- 
eral features which we will discuss below. Each line shown in these plots corresponds to 
a contribution to the total inclusive cross section from a particular combination of quark 
generations indicated in the legends. In fig. HI we show the cross sections involving 

the quark flavor- diagonal a = /3 leptoquark ratios; the vertical dashed lines and horizontal 
arrows indicate the possible location of the r — j- 67 upper bound, as discussed for fig. [31 
Plots for quark fiavor-off-diagonal contributions to the inclusive cross section are shown in 
fig. El 

From these results for the cross sections, we can make several remarks which are generally 
true for most of the scalar leptoquarks. 

1. It is clear that the present limits (2 = 1) on the ratios Xia^^p/Mlq involving first 
generation quarks are more stringent than those limits involving second and third 
generation quarks (these stronger limits are in many cases from the rare process limits 
other than r — ?■ 67 rather than the direct searches at HERA); hence, the allowed cross 
sections for initial state first generation quarks are suppressed despite their larger 
p.d.f. contribution. As a result, with 10 fh~^ and y/s = 90 GeV, it is unlikely that 
the EIC would observe first generation quark e — )■ r events if the leptoquark ratios 
are smaller than half their current limit (the cross section contributions from the 
(a,/3) = (1,1) combination cross the 0.1 fb sensitivity threshold near z = 0.5). To 
probe the first generation quark ratios to the level of one order of magnitude smaller, 
100 to 1000 fb~^ are likely necessary. (However, the cross sections mixing the first 
and second generations in fig. [5] seem to be an exception to this observation.) 

2. For those leptoquark ratios XiaXsp/M^Q involving second and third generation quarks, 
the current limits are less stringent and so the EIC could potentially achieve an order 
of magnitude improvement, or more, with the given energy and luminosity; this is 
possible despite the suppression of these cross sections by the proton p.d.f.s. 

3. We can now address how e — )■ r searches at the EIC would fare in light of the most 
recent r — )■ 67 limit. As fig. H] shows, the present upper limit from HERA and/or other 
rare processes is roughly on the same order of magnitude as the limit on AiqAsq/M^q 
from T ^ e'y for the first generation quarks (the r — )■ 67 upper bounds are near z = 1). 
A factor of two decrease in these limits attained by the EIC would achieve parity with 
the strongest r — )■ 67 upper bound from BaBar. For the second and third generation 
quarks, the r — )■ 67 limits are presently much stronger than the HERA limits (where 
such limits exist — recall that HERA did not set limits for leptoquarks coupling only 
to top quarks in the third generation, and the EIC would likewise be unable to probe 
such couplings). Even the weaker of the two r — 67 upper bounds we calculated seem 
to be generally out of reach of the EIC given our choice of integrated luminosity and 
center-of-mass energy. We emphasize, nonetheless, that these statements apply only 
to the quark flavor-diagonal combinations of scalar leptoquark couplings. Compared 
to present bounds, the EIC sensitivity to quark flavor off-diagonal combinations could 
substantially exceed that of other searches performed to date. 
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FIG. 6. Analogous to fig. SI these plots show the quark flavor-diagonal cross sections for two vector 
leptoquarks. We do not show limits from r — )■ 67 as discussed in the text. 



4. The analysis presented here is predicated on a choice of 90 GeV center-of-mass energy 
and 10 fb~^ integrated luminosity. Larger integrated luminosities will allow the EIC to 
be sensitive to smaller cross sections, and larger center-of-mass energies will increase 
all of the cross sections uniformly since a oc s as in eq. (13. ip . The results of our 
analysis above are broadly improved by such changes. For example, we observe from 
fig. mthat 1000 fb~^ of integrated luminosity would allow the EIC to improve all of 
the leptoquark limits by an order of magnitude or more — even those limits involving 
the first generation quarks, as mentioned above — and surpass the stronger r — 67 
limits for all three generations. 

Although the vector leptoquarks are problematic when calculating the r — > 67 limits, 
the vector leptoquark contributions to e — )■ r are more straightforward. The second term 
of the vector propagator in eq. (14. 5 p could be present in the amplitude for the e — )■ r 
scattering diagrams in fig. [U however, the contribution to the cross sections from this term 
will be suppressed by additional powers of the large leptoquark mass and may safely be 
neglected. Making use of eq. (13. ip . we show in fig. [6] and fig. [7] the quark flavor-diagonal 
and quark fiavor-off-diagonal cross sections, respectively, as a function of z for two of the 
vector leptoquarks, V-^ and Many of the same features discussed above for the scalar 
leptoquarks also generally hold true for the vector leptoquarks (with the principal exception 
that there is no straightforward application of the r — )■ 67 bound to the vector leptoquarks). 
For example, as for the scalar leptoquarks, it is evident that the EIC could obtain an order 
of magnitude improvement on the limits for the second generation quark flavor- diagonal 
ratios AiqAsq/M^q with y/s = 90 GeV and 10 fb~^ integrated luminosity. Finally, we 
point out that the leptoquark V-^ displays a behavior not seen for the other leptoquarks 
presented so far: the cross sections involving_couplings to first and second generation quarks, 
(a, /3) = (1, 2) , (2, 1), are highly suppressedjfl We conclude that the potential of the EIC to 
probe leptoquark-induced LFV(1,3) is mostly independent of whether the leptoquarks are 
scalar or vector particles. 



VI. A MODEL ILLUSTRATION: S^,^ AND NEUTRINOS 



One of the BRW leptoquarks is of particular interest because of its connection to the 
work in 1^. The leptoquark studied by those authors ($& in their notation) has the same 



The source of this suppression is the stringent limit on the decay K 
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FIG. 7. Analogous to fig. EJ these plots show the quark flavor-off-diagonal cross sections for two 
vector lepto quarks. 



spin and gauge group quantum numbers as the leptoquark 3^2- The authors of [l9[ consider 
an SU (5) GUT model in which this leptoquark arises; as a result of the SU{5) symmetry, 
the coupling of this leptoquark to quarks and leptons can be constrained by the neutrino 
sector: by a particular choice of rotation matrices (in particular to avoid proton decay) and 
neglecting phases, the neutrino mass matrix determines the leptoquark couplings. Given 
the experimental constraints on the squared neutrino mass differences and mixing angles, 
the authors perform a scan to determine allowed values of the leptoquark couplings as a 
function of the lightest neutrino mass in both normal and inverted neutrino mass hierarchies. 
The authors introduce a dimensionful coupling parameter Y^-' = x where Fi is the 
leptoquark coupling, equivalent to our A, and Vj\ is the vacuum expectation value of the scalar 
triplet field A which implements the type-II seesaw mechanism. The superscripts i and j 
refer to the quark and lepton generations, respectively (this is a reversal of our notation, Ai^ 
and As^, where the first index refers to the lepton and the second index refers to the quark). 
Their results are shown in figures 1, 2, and 3 of [l^. The quasi-degenerate neutrino mass 
region, where the lightest neutrino mass is O {10^^ eV), is of particular interest because 
of the several orders of magnitude separation in some of the leptoquark couplings Yi. In 
the quasi-degenerate region, satisfying the neutrino mass and mixing constraints forces the 
couplings Y^\ Y^^, and F/^ to be O (10^^ eV), while Y^^ and Y^'^ are O {10'^ eV) and the 
remaining couplings are relatively unconstrained and can take on values between 10"'^ and 
10-1 eV. 

Are there other constraints on this model's leptoquark couplings in addition to those 
imposed by the neutrino masses? As eq. flAll) shows, this leptoquark's coupling-over-mass 
ratios are not constrained by the limit from r — )■ 67 (more precisely, the leading contribution 
to the r — )■ 67 branching fraction from S^i^ is suppressed by m'^/Mf^Q rather than 1/M|q, so 
the upper bounds on AiqAsq/MIq from r — )■ 67 are very weak). Limits from /i — 67 are also 
evaded for the same reason. We instead consider the S^^^ leptoquark's contribution to /i — )■ e 
conversion in nuclei and derive constraints from this tightly constrained LFV(1,2) process. 
The — )■ e conversion rate is calculated from a model-independent effective Lagrangian in 
"sol and can be written in terms of the leptoquark dimensionful couplings Yi as [l9i] 



The quantities V^^^ and V^"'^ are overlap integrals of the muon and electron wave functions 



4 



+2\/("))' . (6.1) 
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with the nucleon densities, calculated for various nuclei in [30| . Using these overlap integrals 



and muon capture rate for gold nuclei [30| and the SINDRUM II collaboration's result for 



the ratio of the muon capture to conversion rates in gold nuclei [8| 



.Au r[l^- + A jZ, N)^e~ + A (Z, N)] uj,„,, ^ 7 n - 1 Q-^^ 

T[^- + A{Z,N)^v + A{Z-l,N + l)] uj^apt ' 



u,apt = 13.07 X lOS-i 



(6.2) 



we can calculate limits on the leptoquark couplings Yi. In the quasi-degenerate region, 
is already constrained to be 10~^ eV as mentioned above; hence, eq. (16. ip and eq. (16. 2p 
dictate that Y^^ < 2 x 10~^ eV , two orders of magnitude smaller than the maximum allowed 
by the parameter scan in [l9| . 

Using the limits on the leptoquark couplings from both the neutrino sector and the /i — )■ e 
bound, we can compute the e — r cross section from eq. (13.11) . replacing Ai^As^/Mlg with 

(M^^q^a)^. After imposing the /i — e limit, the only coupling which is not well- 
constrained in the quasi-degenerate region is Y^^ = Y^"^, for which we take the largest allowed 
value of ~ 10"^ eV. If we take MlqVa = 800 GeV eV as in {va is O (2 el/)), then we 
find the total inclusive cross section for e — >■ r is 0.08 fb for 90 GeV center-of-mass energy 
(well below the total inclusive cross section for this leptoquark calculated using the HERA 
limits). With an integrated luminosity in excess of 10 fb~^, the EIC could be sensitive to 
e — > r events mediated by the S^^^ leptoquark. Therefore, this SU (5) GUT is one interesting 
example of a model which can satisfy existing LFV(1,2) limits on both yU — )■ e conversion 
and — >■ 67 (trivially, for the latter) but still give rise to a LFV(1,3) signal that could be 
observable in the next generation of experiments, specifically e — r conversion at the EIC. 

We conclude this section by noting that closer inspection of the e — r cross section 
reveals an interesting search prospect for this SU{5) GUT leptoquark. The S^^^ leptoquark 
couples to down-type quarks, so the partonic level cross sections featuring a down quark 
in the initial state will dominate due to it from the proton p.d.f.s relative to other initial 
state quarks or antiquarks. An initial-state down quark is only possible via the w-channel 
diagram for F = leptoquarks in fig. [T| hence, the sub-processes proportional to Yl^Y^^ 
will give the largest contributions to the inclusive cross section. Given the neutrino and 
IX ^ e constraints discussed above, this product of the leptoquark couplings is maximized 
for i = 3, corresponding to the leptoquark coupling to an outgoing bottom quark. Thus, 
in e~p collisions, an outgoing r~ and a bottom quark jet may be a unique experimental 
signature of this SU{b) GUT leptoquark. 



VII. CONCLUSION 



We have calculated cross sections for e — >■ r processes as a function of the leptoquark 
ratios AiqAs/j/MIq for both scalar and vector leptoquarks; we have also calculated limits 
on these ratios, where applicable, from the most recent upper bound on r — )■ 67 decay for 
the scalar leptoquarks. Our results, using an EIC center-of-mass energy equal to 90 GeV 
and 10 fh~^ integrated luminosity, indicate that a search for charged lepton flavor violating 
e — J- r events is viable if such processes involve very massive (M^q S> \fs) leptoquarks 
at tree level; furthermore, new limits from the EIC in the event of a null result would be 
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competitive with some of the present r — )■ 67 hmits on quark flavor- diagonal combinations 
of scalar leptoquark couplings, though larger integrated luminosities and energies would be 
necessary to surpass all of the upper bounds imposed by r — )■ 67. On the other hand, an EIC 
search could substantially extend the reach into the space of quark off-diagonal couplings 
compared to current HERA and rare process limits. 

Additionally, we note that leptoquarks participate in strong interactions, and their pair 
production via gluon or quark/antiquark fusion can have large cross sections. Searches for 
leptoquarks at the Tevatron and LHC have been considered [l^, 23, 31-34 1. A search for 



LFV(1,3) leptoquarks at the EIC would be complementary to these leptoquark searches 
at the Tevatron and LHC: the Tevatron and LHC are sensitive to the leptoquark mass 
and branching fractions, whereas e — > r searches at the EIC depend on the ratio of the 
leptoquark couplings and mass, AiqAs/j/MIq. Were the Tevatron or LHC to measure the 
mass of a leptoquark, the EIC could then provide information on its lepton flavor violating 
couplings to quarks and leptons. 

While leptoquarks provide a framework in which e — ?■ r searches at the EIC could be 
feasible, there are many practical questions that are relevant for undertaking this LFV(1,3) 
search. For example, what processes would contribute to the backgrounds? What would be 
the r detection efficiency for the detectors? These questions were previously answered for the 
HERA analyses, and their answers may provide some guidance although the EIC would be 
operating at a different energy and with different detectors. Furthermore, it is important to 
consider the limits on LFV(1,3) from rare processes. From the study of r — > 67 decay in this 
analysis, it is clear that improvements in limits on other rare processes like r — )■ 3e, r — ne, 
and decays of K and B mesons can impact the relevance of e — )■ r searches at the EIC. 
Although 10 fb~^ integrated luminosity could be collected in a relatively short time at 
a high luminosity machine like the EIC, the time at which data collection could begin is 
uncertain and far in the future. Thus, the status of current and future experiments searching 
for LFV(1,3) is necessary information when evaluating the prospective impact of a similar 
search at the EIC. 

Our analysis shows that leptoquark-induced e — )■ r conversion at the EIC is worth further 
study, and the discussion of the SU (5) GUT leptoquark in section IVII shows that searches 
for LFV(1,3) can still be relevant in light of the stronger limits on LFV(1,2) processes. 
Nonetheless, the leptoquark analysis presented here is only an initial study of LFV(1,3) at 
the EIC. There are important theoretical questions to be answered, such as, what are the 
e — >■ r cross sections in other non-leptoquark models with LFV(1,3), and, can searches for 
e — )■ r discriminate between various models? These will be the subject of our future work. 
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Appendix A: Scalar Leptoquark r — )• 67 Results 



Li R 

Here we show complete results for the coefficients contributing to the r — ?• 67 

branching ratio. All scalar leptoquarks are shown. We have neglected the term proportional 
to m^/Mlg in eq. fl4.3p . The charges of leptoquarks in an SU{2) multiplet are differentiated 
using Qo, and 
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